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ABSTRACT 

The effect of gravitational microlensing on the determination of extragalactic distances 
using surface brightness fluctuations (SBF) technique is considered and a method to 
calculate SBF amplitudes in presence of microlensing is presented. With a simple 
approximation for the magnification power spectrum at low optical depth the correc- 
tion to the SBF-based luminosity distance is calculated. The results suggest the effect 
can be safely neglected at present but may become important for SBF-based Hubble 
diagrams at luminosity distances of about 1 Gpc and beyond. 

Key words: gravitational lensing - methods: analytical - stars: statistics - galaxies: 
distances and redshifts - galaxies: stellar content - dark matter 



1 INTRODUCTION 

For many centuries, measuring distances to celestial objects 
has been a central part of astronomy. Historically, building 
, on the trigonometric parallax technique, astronomers have 
' developed a long cosmic distance ladder, within which ev- 
[ ery step is supplied by observing more and more luminous 
r objects further and further away, whose calibration is based 
on distances obtained at preceding steps. 

About two decades ago a new method to determine ex- 
■ tragalactic distances usin g surface brightness fluctu ations 
(SBFs) was proposed bv IXonrv fc Schneider! Jl988l) . The 
method relies on a simple but powerful idea that less nu- 
merous stellar populations show greater relative fluctuations 
in their total brightness as Poissonian fluctuations in the 
number of stars are more pronounced in them compared to 
more numerous ones. This allows one to estimate the true 
average number of stars in and therefore intrinsic luminosi- 
ties of pixels in images of distant galaxies. Comparing these 
estimates with the observed flux then estimates the lumi- 
nosity distance directly. Although the idea of the method 
was clear to some astro nomers long ago (see comments in 
iTonrv fc SchneideJlQsi) . it was not until the mass introduc- 
tion of linear panoramic detectors (CCDs) into astronomical 
observations, that the method could be usefully employed. 

Since then, SBF technique has proven extremely 
fruitful in independent distance determination for galax- 
ies and also in studi es of their stellar popula t ions (for 
review material, see ' Blakcsle e, Aihar fc Tonrvl Il999l or 
iBlakeslee. Vazdekis fc Aihari200lll . Various wavs for the cal- 
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ibration of surface brightness fluctuations - both empirically 
and on the basis of stellar evolution simulations - have been 
developed, allowing to reduce uncertainty in distance deter- 
mination to the level of five to ten per cent, making it one of 
the most precise methods in extragalactic astronomy. The 
range of the method has also increased nearly ten-fold com- 
pared to original estimates, and the distances to galaxies at 
nearly 200 Mpc have be en successfully me asured using the 
Hubble Space Telescope J Jensen et aljl200lfl . 

Since the major limiting factor in the implementation 
of the method is photon noise, the trend in SBF observa- 
tions has been towards the infrared for at least a decade, 
as one gets more infrared photons at a given flux level and 
the SBF signal is dominated by late-type giant stars. In ad- 
dition, IR surface brightness fluctuations turn ed to be very 
interesting for stellar populations studi es fe.g.. |jensen et all 
120031 iMouhcine. Gonzalez fc Lij20o1 . 

One of the external factors which can potentially con- 
tribute to the fluctuations of stellar magnitudes is gravita- 
tional microlensing by compact objects in the Universe. The 
presence of such population has been infe rred from long- 
term photometrical monitoring of quasars jHawkinslllQOa. 
Il996h . The effect of such population on photometrical and 
spectrosco pical properties of quasars has been c onsidered 
previously (ICanizareslll982t iDalcanton et al.lll^^4^ and use- 
ful observational constraints on the properti es of such pop- 
ulation have been obtained JSchneideilll993) . 

This effect has been proposed as a tool to detect 
dark objects in galaxy clusters by observing temporal vari- 
ation s in the flux of a giv e n pixel ([L ewis. Ib ata fc Wvithd 
I2OOOI: iLewis fc Ibatall200ll: iTuntsov et alj|2004ll . However, 
microlensing also introduces extra scatter in flux variations 
among different pixels and makes the distances inferred via 
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SBF method systematically lower than the true ones. The 
correction due to microlensing increases with redshift dra- 
matically, and the compact objects therefore have a poten- 
tial to profoundly affect SBF-based Hubble diagrams at high 
z. 

In the present paper, we investigate how strong this 
effect is and whether it represents a serious issue for the 
present-day SBF observations. Borrowing s ome technique 
from our previous work on galaxy clusters jTuntsov et alJ 
|2004) . we show how microlensing can be taken into account 
for SBF studies, and then calculate relevant quantities un- 
der some general assumptions about potential microlensing 
population. 

We start in Section 2 with an alternative derivation of 
the SBF amplitude which allows to take microlensing into 
account. Section 3 presents a simple summation approxima- 
tion we use for the magnification spatial power spectrum 
and discusses a possible contribution to microlensing by an 
occasional intervening low surface brightness galaxy on the 
line of sight to the source. Then, in section 4 we first show 
that the interstellar correlations within the pixel can be ne- 
glected and then introduce and calculate the correction to 
SBF distances due to microlensing for simulated stellar pop- 
ulations. Section 5 discusses our results. 



2 SURFACE BRIGHTNESS FLUCTUATIONS 
IN PRESENCE OF MICROLENSING 

The method to estimate extragalactic luminosity distances 
d using surface brightness fluctuations (SBF) aims at mea- 
suring the ratio between the dispersion and average value of 
the flux at the pixels comprising a distant galaxy image 



(F) 47rd2 (L) 



(1) 



where L is the total luminosity of the pixel. For a given 
stellar population of the galaxy, the second fraction on the 
right-hand side of this equation can be calculated using stel- 
lar evolution models. Alternatively, one can assume that this 
quantity does not vary much from one galajcy to another 
once they have similar stellar populations in terms of age 
and metalicity and determine its value by calibrating it on 
a sample of nearby galaxies for which the distance D can be 
obtained via other methods. 

Any particular implementation of the method needs to 
find a way to estimate the intrinsic value of ap / {F) from an 
observed CCD image of the galaxy getting rid of difi'erent 
sources of noise present in the real image in comparison to an 
ideal one - non-independence of the pixels introduced by a 
finite- width point-spread function of the optical system, con- 
tamination due to globular clusters of the target galaxy, pho- 
ton count fiuctuations and readout noise, sky background 
etc.. Advanced methods to do so have been developed over 
two decades which have passed since the method was pro- 
posed by Tonry and Schneider (1988) and the intrinsic - in 
the sense that they are unaffected by the sources of noise just 
mentioned - fiuctuations in the fiux can now be measured 
with high precision. 

However, the interpretation of the quantity L has al- 
ways been that it is the true luminosity of the stars com- 
prising a pixel, i.e. the sum of the luminosities of all stars 



the pixel contains. Such an interpretation does not take into 
account that the observed flux contributed by a star can in 
fact be affected by microlensing en route to the observer. In 
this paper we investigate whether this is an observationally 
important issue and in which ways microlensing affects dis- 
tance determination. Since ap/{F) can be determined from 
observations and there is such a simple relation Q between 
this ratio, distance d and a|, / (L) , we will focus on the lat- 
ter quantity and the way these apparent pixel luminosity 
fluctuations are affected by microlensing. 

In the method establishing paper of Tonry and Schnei- 
der (1988), al/(L) was computed in a 'filling number' fiuc- 
tuations manner. That is, it was assumed that there is a 
certain number of different luminosity levels h (minuscule 
letter I is used for individual star luminosities while the cap- 
ital L is reserved for the total luminosity of the pixel), and 
the number Nk of stars of a certain type k in the pixel 
is a Poissonian random variable with the mean parameter 
Nk = N<j>k, where iV is the mean number of stars in the 
pixel and ^j, is the luminosity distribution. This method is 
not applicable when microlensing is involved and we will 
be forced to take a less elegant approach, assuming instead 
that luminosity of a given star is a random variable with 
distribution (pk, and the total number of stars in the pixel is 
a Poissonian variable; microlensing magnification value will 
also be introduced as a random function. It is therefore rea- 
sonable to compare the two types of calculations and point 
out assumptions of the original methods which do not hold 
when microlensing is involved. 

The original gives the following value for the total lu- 
minosity of the pixel: 



L = ^ IkNk. 



(2) 



To find its average value one averages the filling numbers: 

\ k / k k 

where (l)^ is the average individual star luminosity deter- 
mined by its distribution (j). The square of this quantity is 

{Lf = ilk? {Nkf + lklk'{Nk){Nk>) (4) 

k k!^k' 

The average of the luminosity squared is, similarly, 
(L^) = (j2 + ''"^'''^"Nk' \ (5) 

\ fe k^k' I 

k kytk' 

and for the dispersion = {L^) — (i)^ one has 

-i-E('^)'((^^')-<^^)') (6) 

fe 

k^k' 

The value in the brackets on the first line of the right-hand 
side is, for Poissonian distribution, af^^ = Nk ~ Ncpk while 
the second line vanishes as {NkNk') = {Nk'){Nk') because 
the numbers of different stars in the pixel are independent - 
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at least, this is a very reasonable approximation if the galaxy 
stellar population is considered homogeneous. As a result, 
one has 



(7) 



Since Poissonian statistics remains valid for arbitrary low 
mean numbers Nk , the transition to continuous distribution 
density (j>{l) requires no effort. 

The independence mentioned above is lost when mi- 
crolensing is introduced - the magnification values /i inside 
the pixel are correlated: knowledge that a certain number Nj 
of stars have magnification value /ij can be used to predict 
the probability that the number of stars with magnification 
Hji is Nji . We therefore have to resort to a different solution. 

In an alternative approach, the values considered ran- 
dom are 

• the total number of stars in the pixel A*', its distribution 
is assumed Poissonian with mean N; 

• the type of each star, characterized by its luminosity 
h and spatial profile ipi{r) (i = 1,A'^); the values of the 
luminosity I will be randomly drawn with the distribution 

independently of each other^; 

• the coordinates of each star within the pixel with 
i — 1,N , these coordinates will be considered independently 
and uniformly distributed in the pixel; 

• the magnification field realization /i(r), its properties 
will be specified as we progress through the calculation. 

When calculating a\ we will have to average with respect 
to all of these quantities, and, where necessary, appropriate 
averages will be denoted by N , 0, r or /i subscripts at the 
angled brackets, respectively. 

Then, the apparent luminosity of a certain realization 
of the pixel is 



JV 

E 



since 2n%l)i (0) = 1 - it is more convenient to put the normal- 
ization of the spatial profile into the value of U. Then, one 
has 

{L) = {f2{h{tJi),)\ =N{^)^{1)^. (12) 

/ N 

The square of the pixel luminosity is 

AT JV 

= ^^lifiiivi) + ^^lili,^^i{ri)^^i,{ri,). (13) 
One first averages this over different fi{r) realizations 

JV JV 

Homogeneity of the underlying magnification map /i(r) im- 
plies that its autocorrelation function depends only on the 
difference of its arguments: 



{/i(ki)/i(k2))M = P^(ki)5(ki -f ka). 



(15) 



where P^i(k) is the magnification spatial power spectrum, 
or 2n the Fourier transform of the autocorrelation function. 
It follows then that the factor in front of a (5-function in 
the Fourier transform of the cross-correlation between the 
magnifications of stars i and i' is 

P,,.(k) = P^(k)47r2v^,(k)Vi,'(-k) (16) 

and, in particular, the diagonal values of the cross- 
correlation matrix are independent of r: 



(Mi(r))p 



d^kP^(k)lV'.(k)l^ = {^*?), 



(17) 



h/iiiri). 



(8) 



Averaging with respect to Vi is performed by convolving 1141 
with the distribution of r;, i = 1, N. Since these coordinates 
are assumed independent and uniformly distributed in the 
pixel 5 of area jj5||, this leads to 



As the magnification can depend on the type of the source 
star - most importantly, its size, - fi was given the index i 
as well. These values can be calculated by convolving tpi{r) 
and the underlying magnification map fi{r) and therefore 
jii Fourier transform is the product of and /i transforms 
(times 2n as we choose Fourier transformation kernel nor- 
malization in symmetric convention such that only the phase 
is reversed in the inverse transform): 



/i,(k) = /i(k)27rV',(k) 



(9) 



In particular, for a statistically homogeneous magnification 
map the average does not depend on r: 



{/i(k))^ = {M)27r<5(k), 



(10) 



where (/i) is the average maginfication. Therefore, 

{m.)m = (M)M27r^.(0) = (m)m- (11) 



^ For simplicity, we assume a unique relation between li and ipi 
although index i can denote both without any changes to the 
notation. 



where the average value of the correlation function 



{liiti,:)s = \\S\\ J d nd r2 (Mi(ri)Mi'(r2))^ 
can be computed using 11611 : 



d^kp^(k)i/;.(k)7^*,(k)i2^s(k)i= 



(18) 



(19) 



(20) 



Here s(k) is the coordinates characteristic function, or the 
Fourier transform of its distribution density s(r): 



.(r) = ||5| 



res 



(21) 



Since each star type is chosen independently of all oth- 
ers, one has 



= a')^ / d^kP^(k)Kv^(k))^i 

and 



(22) 



4 Tuntsov & Lewis 



{kk' {l^i^^i') s) A 



{1)1 



d^kP^(k)KV(k))2^||2^S(k)l^ (23) 



where the 'luminosity-weighed' averages of the source star 
profile (^'(k))^ and its square {^/^^(k))^ have been introduced 
as, respectively, 



and 



(24) 



(25) 



The number of equal terms in the first and second sums 
of Ijl8|l is A'^ and N{N — 1), and their Poissonian averages 
are N and iV^, respectively. As a result, one has for the 
dispersion of L: 

al = {L^)~{Lf = N{l')^{f^')^+N\l)l (P;^^ - (26) 

with (/J'^)^ and P^'^ defined by the integrals on the right- 
hand side of 1221 and 1231 1. respectively. 

We can change the second term of the last equation into 
a more intuitively clear form by noticing that '!^(0) — (27r)~^ 
and 27rs(0) = 1 while Fourier transform of a constant is 
proportional to a 5-function. As a result, the value in the 
brackets can be calculated as 



piiS 



d2kPAM(k)l{^(k))2||2^s(k)l^(27) 



where Pai^O^) = Pfi(k) — 47r^(^)^5(k) is the power spectrum 
of the magnification fiuctuations A^(r) — /i(r) — (/i)^. 

The dispersion and average value of the observed flux 
in the pixel F differ from those of the luminosity by factors 
{A-Kd^)~^ and (47rd'^)~^, respectively. In the absence of mi- 
crolensing - (fj,) = 1, (fi^) = 1, -P^jf = - this gives the 
following estimate for the luminosity distance d: 

By comparing this formula with 11211 and 1261 one finds that 
when lensing is important this estimate is biased compared 
to the true distance d. The relationship between the two 
values reads: 



■ + A-n- 



Since ^ 1 and P^^ ^ 0, the estimate d is systematically 
lower than the true value d - microlensing introduces extra 
scatter in the pixel fluxes; in addition, mass along the line 
of sight makes distant sources generally brighter on average 
(the first term on the right-hand side of = {^) + 

(A/^^)/(/x)), but this magnification effect is less important 
than the additional scatter (see discussion in section 4) . 

The second term in 129I3L)1 is the most troublesome as 
it could potentially change the recipe for (j\/{L) empirical 
determination - the ratio between the local fluctuation am- 
plitude and mean would no longer be constant and the way 
the mean is subtracted would have to be changed. However, 
as we demonstrate below, typical correlation scale of the 



(29) 



(30) 



magnification map is much smaller than the typical distance 
between the stars. Hence, the correlation term is small. 

Further to this point, the calculation presented above 
does not take into account any correlations which can arise 
between different pixels; taking these correlations into ac- 
count would force us to significantly extend this simple study 
by reconsidering the way in which the intrinsic fluctuations 
in the flux are related to the observed CCD measurement 
statistics. However, interpixel correlations are much weaker 
than interstellar ones and can thus be safely neglected. 

Finally, another important thing that should be ad- 
dressed is the convergence of the variance. Here we only 
calculate the dispersion of the magniflcation value distribu- 
tion, which is an overestimate of the most likely variance of 
a flnite sample of Ns stars forming the galaxy. This happens 
because, for a small source, the dispersion is dominated b' 



a long cc fj, ta il of the distribution density ( Peacock 198' 
ISchneidejlToSTl) . which is in fact unlikely to be properly 
populated by the stars. We can estimate the impact of this 
shortcoming by considering a model magniflcation distribu- 
tion density 



p(m) = 



2(^0-1)^(^^1) 
[(m- 1)2 + (mo -1)2] = 



(31) 



that takes into account three major theoretically established 
properties of the magniflcation distribution function (cor- 
rect normalizati on and the first mo ment p, and the large 
fi behaviour; see iTuntsov et al]|2004l for further details) for 
Ho — 1 = 2{fi — l)/iT. In order to model the finite sample 
effect, we, when calculating the dispersion, will truncate the 
integration not at the value /imax = 2r_B/P, determined by 
the source size R relative to the projection of the Einstein 
radius of the lens 137L but at the value where the cumula- 
tive probability P(> /ijv^ ) drops below the inverse size of the 
sample N~'''. This truncation occurs inside the p(/i) oc 
tail and does not noticeably affect either the normalization 
or the first moment of the distribution. 

In the case when the p.— l ^ 1, which corresponds to the 
low optical depth (estimates in the next two sections show 
that this is the applicable limit), the cumulative probability 
of the magnification greater than /i is 



P(> m) 



l-f- 



2 /i - 1 



(32) 



and therefore for ^jv^ where P(> hn^ ) drops below A'^^ ^ one 
has 



Ativ, =l + -{fi-lWNs-l. (33) 

At redshift z ~ 0.1, the average magnification value due to 
a population of lenses that comprise a fraction fid of the 
critical density, is roughly — 1 ~ 10~'^f2(i (see figures|5||^. 
Therefore, for a target galaxy with A^^ ~ 10^^ stars in it, the 
bound of due to the sample size is fiN^ ~ (10= — 10"^) fid. 
At the same time, the other bound, which is due to the finite 



size of the source is, for a Solar-mass lens, fi„ 



10^-10^ 



depending on the size of the source star. Given that the 
dispersion is roughly proportional to the logarithm of the 
truncation value /i (since p{fj,)fi^ oc fi~^ at large fj,), the 
finite sample effect can decrease the value of the microlensing 
correction to SBF-based distances by a factor of a few. 
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This reduction of the effect magnitude is dependent on 
the number of stars in the galaxy and will not be calculated 
beyond this estimate. As we will see in the following, the 
effect can be safely neglected in present but the finite-sample 
reduction should be kept in mind when the point is reached 
observationally for the effect to be important; in such cases 
the finite-sample reduction should be calculated individually 
for every galaxy. 



of a lens of mass M at the angular diameter distance Dl 
from the observers and Dls from the source is 



AGM DlDls 



Ds 



(37) 



For a cosmological population of compact lenses with 
constant closure density fraction in a fiat Universe , one 
has, at low redshifts z < 0.1 (e.g. iPress fc Gunnlll973l) : 



3 MAGNIFICATION POWER SPECTRUM 

3.1 The summation approximation 

In this section we calculate the magnification power spec- 
trum (k) using an approximation that magnifications due 
to individual lenses simply add up in the limit of very low 
optical depth. 

This is dif ferent from the widely used mu l tiplication 
approximation ijVietri fc OstrikeilllOsi IPeilllOQ,?! : ISchneideil 
[l993) where the magnification /x of a source seen through a 
number of lensing planes is taken as the product of the mag- 
nification values Hi each lensing plane has at the position of 
the source r.: 



(34) 



where the product is over all lensing plane along the line of 
sight. 

We, instead, assume that this value can be calculated 



as a sum: 



(35) 



where the summation extends over all lenses. It is much 
simpler in this approximation to calculate the average and 
the power spectrum of the magnification due to the linearity 
of taking the average. 

This assumption could be justified by saying that the 
two approaches give very similar results when all but one of 
the magnifications are close to unity and a reference to the 
authority of the multiplication approximation. We, however, 
feel the need to go into a little more detail to explain the 
physics behind this assumption. 

In order to do so, one can consider the morphology of 
a macroimage of which a classic example is given in a figure 
by Paczynski (1986). Generally, when there are A*' (point- 
like) microlenses in the field, one has at least A + 1 im- 
ages of which one - the primary - is located close to the 
unperturbed direction to the source and the other A sec- 
ondaries are formed near the lenses. These latter microim- 
ages are highly demagnified unless the line of sight to the 
source passes close to the Einstein-Chwolson ring of one of 
the lenses in which case the magnification gets larger than 
3/\/5 ~ 1.34. The fluxes from all microimages add up. 

When lenses are randomly distributed in space, the av- 
erage number of lenses within their Einstein radii of the line 
of sight is called the optical depth: 



(36) 



r = / AD-KrE{D)n{D), 



AD 



ATvGMn D{Ds - D) 



(38) 



where integration extends from the observer to the source, 
n{D) is the number density of lenses and the Einstein radius 



and thus is relatively small. 

This quantity is also a measure of how many microim- 
ages will be appreciably magnified since the Poissonian 
statistics implies that this number is re~^ which is approx- 
imately r when r ^ 1; we consider just this case of low 
optical depth in the present study. A situation where more 
than one lens contributes noticeable fiux to the macroimage 
is less likely by a factor of approximately 3t/2. 

In addition to always present A -I- 1 image, additional 
pairs of (tertiary?) images may form when the source moves 
into a region bound by a caustic curve. Magnification di- 
verges at these curves so these additional images could con- 
tribute significantly to the total fiux. However, Wambsganss, 
Witt and Schneider (1992) derived the mean number of these 
extra pairs to be Ap ~ in the case of zero external shear 
which is also very low. Granot, Schechter and Wambsganss 
(2003) calculate this quantity for non-zero shear numeri- 
cally and find similar behaviour at low values of r; they also 
show that the statistics of Ap is approximately Poissonian, 
so that T <^ 1 ensures these images can be neglected on 
average. These studies relate to two-dimensional distribu- 
tion of lenses and there have not been studies exploring a 
3D case but at low optical depth there should not be much 
difference. 

Each of the secondary microimages contributes oc [^i — 
l)/2 to the total flux and shears the beam of the primary 
image in such a way that its area, on average, shrinks by 
a fraction (/ii — l)/2. As a result, the input of every lens 
to the magnification value is approximately in — 1 and the 
total magnification is given by (1351 . 

Certainly, ^ and ^ agree in the limit of average in- 
dividual magnification values close to unity. It is also trivial 
that 13411 can be written as a sum simply by applying loga- 
rithm to both sides but averaging this quantity over possible 
magnification map configurations would produce statistical 
properties of In ^ which are hard to convert to those of 
unless some additional assumptions are invoked. 



3.2 Summation for randomly distributed lenses 

We calculate the average magnification and the magnifica- 
tion power spectrum in the approximation 13511 assuming 
that lenses are distributed uniformly and independently of 
each other. 

Consider a cylinder of cross-section C between the ob- 
server at D = and the source at D = Da with A lenses 
inside it. For a point-like lens of mass M at position r; G C, 
Di, individual magnification value is given by (e.g., Schnei- 
der, Ehlers and Falco 1992): 



6 Tuntsov & Lewis 



where 



1 / Ds 



(39) 



(40) 



and the Einstein radius projection onto the source plane is 



-rsi = 



AG Mi D,{Ds-Di) 



(41) 



A V c2 A ' 

where we have assumed that the distance between the i-th 
lens and the source is Dls ~ Ds — Di, which is a good ap- 
proximation for source at low redshift 2; < 0.1 in a spatially 
flat Universe. 

Then for the Fourier transform of = /i — 1 is given by 
the sum 



z>(k) = ^S,(k) 

i 

and each term amounts to 



R2 

= exp 



2tt 



(42) 



(43) 



where we changed to the integration variable — > 
according to l|40|l and introduced a dimensionless Fourier 
transform 



Mi) 



exp[-i$ ■ y] I + 2 _ ^ 



2-K 



R2 



(44) 



r + 2 



- 1 = MO- 



This function is « f ^ for large ^ and tends to unity as 
^ ^ 0; a good approximation can be provided by 

1 



MO 



(45) 



which only slightly (by less than two per cent compared to 
the numerically evaluated C'oiO) overestimates its value at 
low ^ < 1. 

Averaging with respect to the position of lens in the 
appropriate plane G £ produces, in the limit of infinite 
cylinder cross-section C —* Ti? 



exp 



\c\ 



and if the distribution of lenses along the D coordinate is 
some lpd{D) oc n{D), one has in v(k) a sum of TV equal 
terms: 



27r(5(k) 



|£|| / dDn(D) 



dD 7rr|(7:))7i(D)i>o(ri3k).(47) 



Since i>o(0) = 1 and the average number of lenses within the 
cylinder is iV = J dD \ \C\\n{D) we arrive at 

{i>(k))f,. = {,y)^2TvS{k) = ((m)m - 1) 27r5(k) = 2r-27r5(k).(48) 



The optical depth is proportional to the first power of mass 
of the lens (r_E oc \/M), therefore a distribution of masses 
can be taken into account by using the average value of the 
mass. 

In a similar fashion, when calculating the correlation 
function, one has: 

JV N 
f>(ki)^(k2) =^^.(ki)^,(k2)+^!>,(ki)/>,,(k2) (49) 

N 

= ^ r%i exp [-i (ki + k2) fi] vq (r^iki) vo (f Eik2) 

i 
N 

+ ^ f%if%i, exp [-ikif i - ik2f j'] uq (rEi^i) i>o {fEi''i^2) , 

where for typographical clarity we used lens coordinate pro- 
jections onto the source plane: fi = ViDs/Di. 

Averaging of the second sum produces N{N — 1) equal 
terms of the form: 

(i>,(k))^ = (f ) ' 4^'5(ki)5(ki -f k2) (50) 
while the first sum produces A'' terms of the form 

(z>,(ki)/>,(k2))M = (51) 

27r / D.^ ^ 



5(ki + k2) ( ^ (gj) r|,(A) 1^0 (kifB.(A))r 



Taking into account equation 148 > and the fact that the 
average values of A'' and A(A'' — 1) are N and 7V^, respec- 
tively, one has 

(i>(ki)^'(k2))p = (5(ki + k2) (PA.(ki) + 47r"5(ki)(i.)^) ,(52) 
where the power spectrum PAi^(k) = PAfj(k) is 

PAM(k) = j AD2-Kn{D)(^^ f%{D)\vo{'kfE{D))\^ (53) 



= 2 



dr(D)f|(D)!i>o(krE(D))i = 



- that is, twice the cumulative optical depth times the differ- 
ential optical depth- weighed average value of f%\v(,{\y.f e)]"^ . 

At low redshifts, using 145L one can write the follow- 
ing approximation for a population of lenses with constant 
number density: 



Pf,(k) = {l + 4.T)A-K^5{k) + 



9(l + fcfE+2fc2f|) 



+ o(r),(54) 



valid to within four per cent for any k, where the Einstein 
radius projection for the lens halfway to the source is 



Te 



AGM 



(55) 



3.3 The contribution of LSB galaxies 

An inevitable source of noise to the microlensing correction 
that needs to be quantified is the potential presence of a low 
surface brightness (LSB) galaxy on the line of sight to the 
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target galaxy. Even if the dark matter were not composed of 
any compact sources, such a galaxy would cause microlens- 
ing on its stars and as up to a half of t he galaxy population 
is in the LSB form (Bothun. Impev fc McGaugh 1997). a 
significant fraction (~ 10 — 50 per cent) of the potential tar- 
gets at high redshifts z > 1 will be lensed by these galaxies, 
although a confident estimate is difficult to make because 
of our still poor knowledge of the density and clustering of 
galaxies at the low end o f the surface brightness distribution 
Jimpev fc BothuJl997l) . 

The effect of such a galaxy on the microlensing contri- 
bution to the SBF amplitud e can e asily be estimated using 
the results of iTuntsov et all i2004r) . where some of the for- 
malism of the present paper has already been presented, and 
can be straightforwardly borrowed to consider the effect of 
a two-dimensional lens on the fluctuation properties. Using 
their equation (14) and neglecting the interstellar correla- 
tion term in l|26^ of the present paper (see below) we find 
that the effect of a low surface brightness galaxy on the SBF 
amplitude can be expressed as an extra factor of 



(56) 



where ef, is the 'cor relation amplitude' defined by (19) of 
iTimtsov et all t004\ . 

To determine the two factors above, we need to know 
the typical mass of the microlenses and calculate the di- 
mensionless convergence k and shear 7 of their population 
in th e foreground galaxy, as well as their 'effective' values 
jPacz vrlski 1986; Kavsor, Rofsdal & StabcU 1986) given by 
the mix of the compact and smoothly distributed matter in 
the galaxy. As for the first ingredient, the resulting value of 
the dispersion is very weakly dependent on the ratio of the 
microlens mass provided that the ratio of its Einstein radius 
to the source size i s large enough (see the Appendix B of 
iTuntsov et al.ll2004l and discussion following (1671 in this pa- 
per), which is the case here, as shown in the next subsection. 
We, therefore, can adopt solar values for an estimate. 

Then, one has — |(1 — k)^ — 'y^\~^ and 20 

e?,(K.ff.7.ff,f fl/jt;) calculated in subsection 2.2 and shown ^ = lOOpc^'ioo 60.2 = 3.1 x 10 cm7:iioo fo.; 



a dark and smooth form dlmpev fc BothuJl997fl . the total 
convergence is Ktot ~ 10~^ — 10~^; the difference between 
the stellar and effective convergence values is insignificant 
in this case. 

To evaluate the shear 7, one would need an accurate 
determination of the surface brightness density over the en- 
tire surface of the foreground galaxy, which is not possible 
if it is invisible. Generally, however, the scaled shear takes 
the values of the same order as the scaled convergence, and 
~ 0{k) in this case (again, assuming the typical Ein- 
stein radius is much larger than the typical target star size). 
Therefore, a typical contribution to the correction of the 
SBF amplitude due to an occasional intervening LSB galaxy 
along the line of sight is of order 10~^ — 10~^ at moderate 
redshifts. This is less than what we find in the next section 
for a cosmological distribution of compact objects if they 
make up any significant fraction of the matter content of 
the Universe, but can be comparable and even dominant in 
case no such distribution is present. 



4 EVALUATION OF CORRECTIONS 

4.1 Unimportance of interstellar correlations 

The expressions above are sufficient for the calculation of 
{^J.^),/, and -P^jf once the average source star profiles {?/>(k))0, 
and the pixel shape s(k) have been specified. 
However, we will first make some estimates to see what 
the typical values of the parameters of the problem are. 
There are four main length scales of the problem: the source 
size R, the Einstein radius projection to the source ve and 
the pixel size x, as well as the typical projected separation 
of stars in the galaxy A. 

The pixel size is determined by the resolution of the 
CCD matrix and the distance to the source. The typical 
pixel size is, in the angular measure, of order 9 « 0.2" « 
lO"'^ rad. At the distance D this translates into: 



(58) 



in Figure 1 of lTuntsov et al.l ll2004) . The value of the scaled 
convergence can be calculated from the observed surface 
brightness of the foreground galaxy u (in magnitudes per 
square arcsecond) and the mass-to-light ratio T of the com- 
ponent of interest by comparing it with the bolometric irra- 
diance of the Earth by the Sun as 



= 4.1 kpc 2 6*0.2 = 1.3 X 10 cm 2 610.2, 



^ttGMq DlsDol /Irad^ 



AU" 



Dos 



(^) ( 



l + 20''l0°-*('"o'-''{57) 



2.04 X io-3-o-*{"-22-^)t(1 + zA^ 



Dos 100 Mpc ' 

Here Dls, Dql and Dqs are the angular diameter dis- 
tances between the lens and the source, the observer and 
the lens and the observer and the source, respectively, and 
zi is the lens redshift. Given that the mass-to- light ratio 
for stellar populations are of order unity (this value is not 
known well for LSB galaxies though), and for a galaxy to 
be non-detectable one needs its surface brightness to be be- 
low ~ 23™ — 24"* per square arcsecond, the typical values 
of the scaled convergence in stars are of order 10""^ — 10~^ 
up to moderate redshifts. If, in addition, one assumes that 
a significant amount of mass in LSB galaxies is present in 



where Dioo = D/lOO Mpc, 6I0.2 = 61/0.2" and the second line 
is valid at low redshifts z < 0.1. 

The Einstein radius projection (for the lens halfway to 
the source) is, according to II55|I . 

T-E = 4.5 X 10"^ PcD\^l rn''^ = 1.4 x 10^^ cmD^^o rn^\^9i) 

no in-2 1/2 1/2 01 inl6 1/2 1/2 

= 2.8 X 10 pcz ' m ' = 8.7 x 10 cmz ' m ' , 

where m = M/Mq. 

The stars which contribute non- negligibly to the lu- 
minosity all lie on the main sequence or above it in the 
Hertzsprung-Russell diagram. Their sizes therefore lie in the 
range: 

i? 3 X 10"^ pc - 3 X lO"** pc « 10^" cm - lo" cm. (60) 

The projected density of stars in galaxies hardly exceeds 
10* pc~^ in their densest parts and is normally of order 
lO'^ pc~^ on average. Therefore, for the mean separation of 
stars in the pixel one has 

A > 10"^ pc - 10"^ pc « 3 X 10^*^ cm - 3 X lo" cm. (61) 
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Therefore, in all extragalactic cases one has R <^ fs ^ ^■ 
The respective scales in the Fourier images domain, which 
are given roughly by the inverse of these quantities and rep- 
resent size of the region in k plane within which ^p, Pa/j and 
s differ appreciably from zero satisfy the opposite relation. 

This observation allows us to interpret and estimate the 
second term in the brackets of H30^ . Namely, the integral l|27|l 
is contributed almost exclusively by the region of size 1/x 
within which the values of Pa^i and V are nearly constant 
and equal to 4rf|;/9 and 1/4%^, respectively. As a result, 
we can write the following estimate for the correlation term 
using H54^ : 



An 



-irf% 
97r2 ■ 

3^5 



(62) 



Therefore, NP^^ is proportional to the average number of 
stars within the Einstein-Chwolson circle of the typical lens 
Ne- Therefore, taking into account the estimates 1591 
ffs ~ 10~ — 1 for distances from 100 Mpc up to a few Gpc. 
Finally, for the luminosity ratio one has: 

1 r,0.4{M-(M>) 



(63) 



where, in accordance with the tradition adopted in SBF 
studies, I has been used to denote the luminosity-weighed 
average luminosity of source stars (this quantity is also 
called the surface brightness fluctuations amplitude in view 
of 12811 '): upright M are used for corresponding absolute mag- 
nitudes to avoid confusion with the lens mass. Due to a great 
range in stellar luminosities, M assumes values close to the 
absolute magnitudes of red giant branch stars while (M) are 
closer to those of red dwarfs dominating the main sequence. 
The ratio of their luminosities ranges from about a hundred 
to thousands, depending on the band. Therefore, 



N 



2 p; 



:1pS 



tNe 



]^q0.4(M-(M)) ^ ^ 



{P) (m) 9n^f,) - - ■ ^^^^ 

The other term in the brackets of ^M)\ is given by the convo- 
lution of P^ with the source star profile. Assuming all stars 
are Lambert discs of the same radius P, one has 



7^(k) = ^(fc) = 



Jl(fcP) 



TxkR 

and 1171 1 taken with 1541 1 gives 
where 



2 + 



97r-^"(p) 



fa{p) = P 



C(1 + C + 2C2) 8 



(65) 



,(66) 



(67) 



the approximation is valid for large p. 

At D = 100 Mpc, ffi/P ^ 2 X (10"* - 10^) for solar- mass 
lenses and source star radius P = (10 — 100) P©, typical 
for red giants which normally dominate SBF signal. Thus, 
the first term in 13UII is much greater than the second for 
all sensible values of parameters. This supports the claim of 
the paragraph following that equation, and we will neglect 
the second term hereafter. 

The only conceivable situation where the correlation 
term is non-negligible is if the microlenses are much more 



massive than the Sun. Since the average number of stars 
within the Einstein ring of the lens is proportional to the 
mass of the star, we see taking into account 16411 and the 
estimate for the luminosity ratio, that a lens mass of order 
(lO"' — 10^) Mq is needed for the correlation term to be- 
come of order r. Increasing the mass would also affect the 
average in 1661 but the dependence is much weaker - nearly 
logarithmic - in this case. 

When the correlations are important the whole method- 
ology of the SBF method breaks down as the dispersion of 
the pixel luminosity is no longer proportional to the aver- 
age. Such change would manifest itself in the statistics of 
the surface brightness fluctuations generally enhancing the 
probability of large deviations from the mean value. Careful 
study of the observed fluctuation statistics therefore presents 
a potential means to detect a population of massive compact 
lenses. However, results of the next subsection show that the 
effect is unlikely to be important and we therefore do not 
extend our study beyond the present observation. 



4.2 Corrections in cosmological settings 

With interstellar correlations seemingly negligible, we can 
focus on a single quantity 



= l + 2r+ / d'kPAM(k)KV''(k))j+5(r). (68) 



In order to evaluate this quantity three inputs are needed - 
the optical depth to microlensing r and magnification fluc- 
tuations power spectrum PAfj(k), given by 1361 and 15311 . 
and the squared-luminosity-weighed average profile (iA^(k)) 
defined by ll^ . 

The former two quantities depend on the global prop- 
erties of the lensing population and include underlying cos- 
mology through the use of angular diameter distances. In 
order to evaluate the m, the approach of a classic paper by 
IPress fc GunnI il973t) will be used. For reasons of the econ- 
omy of hypotheses we also assume that the number density 
of compact lenses is constant in comoving coordinates and 
all lenses have the same mass M while for the evaluation 
of distances the now-standard cosmographical model with 
A = 1 — SI = 0.74 and a flat geometry will be used. 

The angular diameter distance D to redshift z is then 
D{z) = cr{z) / Ho{l+z), where Ho is the present-day Hubble 
constant and 



r{z) 



Az' 



^^{1 + z')'i + 1-0. 



(69) 



If the present-day value of the critical density in compact 
objects is Q,d and it increases into the past oc (1 -I- z)"^, one 
has {Dls = Ds- Dl(1 + zl)/{1 + zs)) 



3 f dz{l + z)r{z) 

2 J Jnn +zy' + i-n 





1 



r{z) 



r{zs) 



while the magniflcation fluctuation power spectrum 
dz{l + zf[r{z,)-r{z)]'' 



PAM(k) = Si^dTH 



(1 + zs)^^/U{rTzfTT~Q 



(70) 



(71) 
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Figure 1. Magnification fluctuation dispersion as a function of 
stellar radius (in solar units) for different redshifts. The values 
are computed for a cosmological model described in the text and 
8.6 X 10^^ cm corresponding to a lens mass of 1 Mq. 



2 2 (l + z)r(z,) 



r{z)] 



where the approximation 14511 has been used and the cos- 
mological Einstein radius rn introduced as 



4GM 
cHo ■ 



(72) 



The integrals above need to be calculated numerically. 

The only ingredient left is the average source star pro- 
file {i/''^(k)). As the contribution of each star is proportional 
to the square of its luminosity, it is even more biased to- 
wards highest-luminosity giant stars than I. However, the 
contribution of smaller stars here cannot be neglected be- 
cause PAfi{k) decreses slowly with k and small stars have 
more extended ip profiles. We therefore need to consider a 
model population of stars and find an average profile in this 
way. 

We will assume that all stars are Lambert discs so that 
their spatial profiles are described by I65II and only differ 
in R, which is independent of the band; we do not consider 
any limb darkening in this simple study. Luminosity, on the 
contrary, depends on the wavelength, and therefore (i/''^(k)) 
differs from one spectral band to another which makes {/i^)^ 
depend on the band, too. Since gravitational lensing itself is 
achromatic, dependence via profile is the only way in which 
spectral bands enter {/^^)</,. 

To compute ('!/)^(k)) we use the Bag of Stellar Tracks 
and Isochrones (BaST I ) population synthesis software 
JPietrinferni et al.1 1200I ICassisi et all l2005l ICordier etall 
l200^) ^In addition to stellar luminosities in nine bands - 
U, V, B, R, I, J, H, K and L - the code computes bolomet- 
ric luminosities and effective temperatures allowing us to 
calculate their radii using Stefan-Boltzmann law. The code 
has been used to compute SBF amplitudes in the infrared 
and wa s shown to agree with 2MASS data for Magellanic 
Clouds jGonzalez. Liu fc Bruzuaill2005^ ■ 



The code is available at \http://www.te. astro.it/BASTI/\ 



For the economy of hypotheses we have chosen to use 
Galactic bulge star formation history to represent that of 
the target galaxy stellar population. A number of simula- 
tions with different initial seed values were performed for 
about a million stars each. They do not show any noticeable 
difference in terms of (n'^) 4, / (fj) from each other. As a com- 
parison reference, the same values for the Galactic disc-like 
populations were computed. 

Since both r and (o"^)^ are proportional to the density 
of the microlenses Qd, we present our results in terms of the 
correction q = (2t -I- (Aii^)^)/Q,d defined such that the true 
luminosity distance to the source galaxy d and its estimate 
d are related by a form of 131)1 : 



d'iz) 



d\z) 



l+q{z)^d 



+ o(r). 



(73) 



Note, that when the microlenses represent a dominant 
fraction of the Universe matter density, the correction de- 
fined above is the one relative to the so-called 'empty beam' 
distances (Dyer & Roeder 1972; 1973). The combined action 
of all the masses in a respective 'filled beam' would produce 
an average magnification of (l-|-4r) and therefore correction 
to the 'filled beam' distance is g — 2r = {Afj?)^. The differ- 
ence between the two values is less than twenty per cent in 
the situations considered below. 

To calculate q, we first compute the dispersion of the 
magnification as a function of stellar radius for a given 
source redshift Zs using approximation H71|l and source pro- 
file 165|l - figure Q shows examples of such functions - and 
then average it using results of numerical simulations sepa- 
rately in each spectral band. Then, the 2r value is added. 

Figure |5] shows the results of our calculations. Cor- 
rection q, defined above, and the optical depth r/Qd are 
shown as functions of the redshift z and true luminos- 
ity distance d (upper and lower horizontal scales, respec- 
tively). We adopted the currently favoured spatially flat 
Universe for this calculat ion with parameters f l = 0.26 and 
Ho = 73km-s"^ -Mpc"^ ijEidelman et al.l2004h . as well as a 
microlens mass of 1 Mq. An old. Galactic bulge-like popula- 
tion has been used to calculate the microlensing correction. 
We reiterate that the finite-sample effect discussed at the 
end of Section 2 should be considered for each galaxy indi- 
vidually, when the effect is computed for applications to real 
observations; it further reduces the magnitude of the lens- 
ing correction to SBF-based distances and depending on the 
number of stars in the target galaxy, this reduction can be 
a factor of a few for smaller galaxies. 

As shown in the previous subsection, typical Einstein 
radii of the lenses are much larger than the stellar radii at 
most of the distances of interest and therefore for lenses of 
much larger mass the relation H6()^ in the large p limit of l|67^ 
holds; at the same time, fs oc y/W. Therefore, higher masses 
would produce corrections that are greater by an additive 
term of roughly Tln{M /Mq)/{18tv), which is just 12 per 
cent of the 2r value for lens masses as large as 10® Mq . 

In order to avoid crowding in the graph, only bolomet- 
ric, ultraviolet (U), visual (V) and two infrared (/ and H) 
values of q are depicted. In fact, B values are very close to 
those for V, being about five per cent lower than the latter 
for 2 < 0.08 and slightly higher (by up to four percent) for 
larger redshifts. K and L graphs would be indistinguishable 
from those of H on the scale of this graph, while R and J 
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Figure 2. Correction q, defined by 1731 and the optical depth 
r/Qii as a function of redshift z and the (true) luminosity distance 
d. Only the bolometric value and four filters - U, V, I and H are 
shown to avoid crowding in the graph. A spatially flat (f2 + A = 1) 
Universe with Q = 0.26 is assumed and Hq = 73 km ■ s^^ • Mpc"^ 
are adopted. 
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Figure 3. Same as figure |^ computed for a Galactic disc-like 
stellar population. The two graphs are very similar. 



take values roughly half-way between V and /, and / and H, 
respectively. 

Figure 13 shows the same quantities obtained for a disc- 
like population. One can see, that unlike the SBF fluctua- 
tion amplitude, the magnification correction value q is very 
weakly dependent on the stellar population. The comment 
above concerning the finite-sample effect can be repeated in 
full here. 

The complex behaviour of the graph in terms of the de- 
pendence on colour at low redshifts can be explained by the 
non- linear form of (A/i^)(7?) dependence. It is proportional 
to the logarithm of te/R for small source radii and falls of 
as for much larger sources, which introduces an extra 

weighing factor in averaging for close galaxies, when a signif- 



icant fraction of the optical depth is contributed by nearby 
lenses. This behaviour is less pronounced for the disc-like 
population of figure |3] 

By comparing the average values of q with (A/i^)(7?) in 
figures|21andQwe see that the typical effective radius of stars 
that contribute most to the correction is of order ten solar 
radii, except short-wavelength bands U, B, V, where one 
can see a significant contribution due to stars with radii of 
solar order. Therefore, as is the case with the surface bright- 
ness fiuctuations amplitude, the correction is dominated by 
late-type giant stars. This is quite natural because the aver- 
age value we are interested in is weighed with the square of 
stellar luminosities, making red giants the most important 
contributors - especially, for old stellar populations such as 
the one used in our calculations. 



5 CONCLUSIONS 

In the present paper we have considered the correction to 
distances inferred from the surface brightness fluctuations 
analysis, caused by a potential population of compact ob- 
jects along the line of sight to the source galaxy. Using a 
simple approximation for microlensing magnification fluctu- 
ations we have studied how this effect can be taken into 
account and whether such accounting is needed in obser- 
vations. Although we have considered microlensing as the 
main factor of individual star brightness fiuctuations, the 
technique developed should be applicable to a wider range 
of geometrical variability. 

By combining the results shown in figure |5| with the 
upper limit of ^ 0.3 on the total fraction of matter in the 
energy budget of the Universe, one can see that microlensing 
can for now be safely neglected when interpreting SBF-based 
distances. The present-day observational facilities limit the 
range of the method to less than 200 Mpc while the uncer- 
tainty of this method is presently between five and ten per 
cent. The finite-sample effect discussed at the end of Section 
2 further reduces the magnitude of the lensing correction. 

However, for the next generations of telescopes, this sit- 
uation can - at least, in principle - change. The main factor 
limiting the accuracy of the surface brightness fluctuations 
method is the photon noise, which is of geometrical origin. 
Therefore, a naive expectation would be that increasing the 
diameter of the telescope n-fold would lead to a directly 
proportional increase in the range in terms of the comov- 
ing distance c/Hor{z). Taken with the estimates given in 
this study, this implies that microlensing might just start 
to matter when space telescope diameters get closer to a 
ten-metre scale. 

Alternatively, one might notice that q turned to be of 
the same order as the total optical depth. Therefore, any 
large overdensity of compact microlenses - such as those 
potentially associated with galaxy clusters - would intro- 
duce noticeable change in the surface brightness fluctuations 
statistics and can thus be detected in this manner. However, 
the properties of strong lensing mapping need to be known 
in full detail for such project since local variations of the 
apparent luminosity distance are greatest near the critical 
curves of the mapping which yield galaxies most suitable for 
SBF studies. Temporal variability, therefore, seems to be a 
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better tool to study compact dark matter in galaxy clusters 
jTuntsov et al.ll2004^ ■ 

We conclude that microlensing by a cosmological pop- 
ulation of compact objects is unlikely to significantly bias 
determination of distances using surface brightness fluctua- 
tions technique. Similarly, present accuracy of this method 
does not allow to put useful constraints on the properties of 
such population from the absence of systematic discrepan- 
cies between distances determined via SBF and other meth- 
ods. The effect considered in this paper, however, should 
not be overlooked when the accuracy and range of the SBF 
method reach the level of a few per cent and a gigaparsec, 
respectively. 
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